Chemical Abundances in the Universe: Connecting First Stars to 
Planets 

Proceedings IAU Symposium No. 265, 2009 © 2009 International Astronomical Union 

K. Cunha, M. Spite & B. Barbuy, eds. DOI: 00.0000/X000000000000000X 

Fe I/Fe II ionization equilibrium 
in cool stars: NLTE versus LTE 

Lyudmila Mashonkina 1 ' 2 , Thomas Gehren 2 , Jianrong Shi 3 , 
Andreas Korn 4 , and Frank Grupp 2 

1 Institute of Astronomy, Russian Academy of Science, 
Pyatnitskaya 48, 119017 Moscow, Russia 
email: lima@inasan.ru 

2 Universitats-Sternwarte Munchen, 
Scheinerstr. 1, 81679 Munchen, Germany 
email: lyuda, gehren, fug@usm.lmu.de 

3 National Astronomical Observatories, Chinese Academy of Sciences, 
A20 Datun Road, Chaoyang District, Beijing 100012, PR China 
email: sjr@bao.ac.cn 

4 Department of Physics and Astronomy, Uppsala University, 
Box 515, 75120 Uppsala, Sweden 
email: andreas.korn@fysast.uu.se 

Abstract. Non-local thermodynamic equilibrium (NLTE) line formation for neutral and singly- 
ionized iron is considered through a range of stellar parameters characteristic of cool stars. A 
comprehensive model atom for Fe I and Fe II is presented. Our NLTE calculations support the 
earlier conclusions that the statistical equilibrium (SE) of Fe I shows an underpopulation of Fe I 
terms. However, the inclusion of the predicted high-excitation levels of Fe I in our model atom 
leads to a substantial decrease in the departures from LTE. As a test and first application of the 
Fe I/II model atom, iron abundances are determined for the Sun and four selected stars with 
well determined stellar parameters and high-quality observed spectra. Within the error bars, 
lines of Fe I and Fe II give consistent abundances for the Sun and two metal-poor stars when 
inelastic collisions with hydrogen atoms are taken into account in the SE calculations. For the 
close-to-solar metallicity stars Procyon and /3 Vir, the difference (Fe II - Fe I) is about 0.1 dex 
independent of the line formation model, either NLTE or LTE. We evaluate the influence of 
departures from LTE on Fe abundance and surface gravity determination for cool stars. 
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1. Introduction 

Iron plays an outstanding role in studies of cool stars thanks to quite numerous lines in 
the visible spectrum, which are easy to detect even in very metal-poor stars. Iron serves 
as a reference element for all astronomical research related to stellar nucleosynthesis and 
chemical evolution of the Galaxy. Iron lines are used to determine the surface gravity, 
log g, and the microturbulence £ of stellar atmospheres. In the atmosphere with T c ff > 
4500 K, neutral iron is a minority species. The ionization equilibrium between Fe I 
and Fe II and the excitation equilibrium of Fe I easily deviate from thermodynamic 
equilibrium. Since the beginning of the 1970s a number of studies attacked the problem 
of non-local thermodynamic equilibrium (NLTE) for Fe (e.g., Athay & Lites (1972), 
Thevenin & Idiart (1999), Gehren et al. (2001)). However, a consensus on the expected 
magnitude of the NLTE effects was not reached. 

In this study, we update the model atom of Fe I-II treated by Gehren et al. (2001) 
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(hereafter Paper I) and apply it to analysis of the Fe spectrum in the Sun and selected 
cool stars with the aim of empirically constraining the role of inelastic collisions with 
hydrogen atoms in the SE of Fe I-II. 

2. The Fe model atom 

In all previous NLTE calculations, the model atom of Fe I was build using measured 
energy levels. The experimental analysis of Nave et al. (1994) with later updates provided 
965 energy levels for Fe I. A comparison with the calculated Fe I atomic structure (Kurucz 
(2007)) reveals that the system of measured levels is nearly complete below excitation 
energy, E exc , 5.6 eV, however, laboratory experiments do not see most of the high- 
excitation levels with E exc > 7.1 eV. As already shown in the first NLTE studies, the 
main NLTE mechanism for Fe I is the overionization of low-excitation levels by ultraviolet 
radiation. The role of high-excitation levels is to compensate, in part, for population losses 
via collisional coupling to the large continuum reservoir, with subsequent spontaneous 
transitions down to low-excitation levels. Therefore, the system of levels in the model 
atom of Fe I has to be fairly complete at least up to 0.5 eV (mean kinetic energy of 
electrons in the atmosphere) below the ionization limit. 

For Fe I, our model atom was constructed using all known energy levels and the 
predicted levels with E exc up to 7.83 eV, in total, 2970 levels. Multiplet fine structure 
was neglected for all terms. The predicted and measured levels with close energies were 
combined resulting in 233 terms. In addition, six super- levels were made up from the 
remaining predicted levels. For 11958 radiative transitions occurring in this atom of 
Fe I, (//-values were taken from the Nave et al. (1994) compilation, where available, and 
Kurucz (2007) calculations. Photoionization cross-sections of the IRON project (Bautista 
(1997)) have been used for 149 levels and a hydrogenic approximation for the remaining 
levels. The collisional rates were computed as in Paper I. 

For Fe II, we rely on the reference model atom treated in Paper I. In this study, it was 
reduced and includes now the levels with E eyic up to 10 eV. The main uncertainty of the 
NLTE calculations for Fe I and II is the treatment of the poorly known inelastic collisions 
with hydrogen atoms. We employ the formula of Steenbock & Holweger (1984) for allowed 
transitions and a simple correlation between hydrogen and electron collisional rates, 
Ch = C e y/ (m e /riiH)NH/N e , for forbidden transitions. Calculations were performed with 
the hydrogen collision enhancement factor Sn, which was varied between and 3. 

3. Results 

The coupled radiative transfer and statistical equilibrium equations are solved with 
an improved version of the DETAIL program (Butler & Giddings (1985)) based on the 
accelerated lambda iteration. All calculations are performed with plane-parallel, homo- 
geneous, LTE, and blanketed model atmospheres computed with the MAFAGS-OS code 
(Grupp et al. (2009)). 

For comparison with observed data, a total of 43 lines of Fe I and 18 lines of Fe II 
were chosen. For the Sun and HD 84937, the analysis was extended to a larger line list 
including 271 lines of Fe I and 34 lines of Fe II. The Sun is also used as a reference 
star for a line-by-line differential analysis of stellar spectra. Solar flux observations were 
taken from the Kitt Peak Solar Atlas (Kurucz et al. (1984)). The absolute solar iron 
abundances were determined using (//-values from O'Brian et al. (1991) and Melendez 
& Barbuy (2009) for Fe I and Fe II, respectively. We find that virtually all models of 
line formation, whether LTE or NLTE with Sh ^ 0.1, lead to acceptable solar ionization 
equilibria within their ler error bars. To show the maximal NLTE effect on abundance 
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Table 1. Stellar parameters and iron abundances obtained for selected stars 
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0.9 
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1.4 


0.04±o.O3 
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determination, Table[T]presents the average abundances for both ionization stages derived 
from the NLTE with S H = (denoted as NLTE ) and LTE calculations. 

Four stars with effective temperature and surface gravity measured from the model- 
independent methods were chosen to investigate the ionization equilibrium between Fe I 
and Fe II for various Sh values. They are listed in Table[T] together with the T e s and 
log g values taken from Di Folco et al. (2004) for HD 10700 (r Cet), Allende Prieto et al. 
(2002) for HD 61421 (Procyon), Korn et al. (2003) for HD 84937, North et al. (2009) for 
HD 102870 (/3 Vir). Observational data were obtained with the FOCES spectrograph at 
the 2.2m telescope of the Calar Alto Observatory during a number of observing runs 
between 1997 and 2005, with a spectral resolution of R ~ 60 000 and a signal-to-noise 
ratio S/N ^ 200. 

The NLTE, Sh = and LTE abundances obtained from the lines of Fe I (denoted 
as [Fe/H]i) and Fe II ([Fe/H]n) are presented in TableQ] It is worth noting that, with 
the updated model atom of Fe I-II, the departures from LTE are substantially smaller 
compared to those from the previous studies. For example, with Sh — 0, we obtain 
an average NLTE abundance correction Anlte = log £nlte — log £lte = 0.22 dex for 
the Fe I lines in HD 84937, while the corresponding value amounts 0.40 dex in Korn 
et al. (2003). FigureQ] displays the abundance difference between Fe I and Fe II for 
various assumptions for the hydrogen collisions. We find that NLTE with pure electronic 
collisions (Sh = 0) is not acceptable for HD 84937 and r Cet. This indicates the need 




Figure 1. The difference in abundance between Fe I and Fe II in selected stars from calculations 
with various line formation models. For each star, the error bars is indicated in the upper part 
of panel. 
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for thermalizing processes not involving electrons in the atmosphere of metal-poor stars. 
For each object, the NLTE effect on abundance determination is small (within the error 
bars) when hydrogen collisions are included with Sh 7? 1- For Procyon and Vir, the 
mean Fe abundance from Fe 1 lines is about 0.1 dex lower compared to that from Fe II 
lines. The origin of such a discrepancy will be investigated in a forthcoming paper. For 
the Fe I/Fe II ionization equilibrium in two metal-poor stars, LTE seems to be as good 
as NLTE with S H > 1. 

The NLTE calculations were performed with Sh = 1 for the small grid of model 
atmospheres with T c s = 5000 and 6000 K, [M/H] = —1 and —3, and log g ranging 
between 2 and 4 in order to inspect the departures from LTE depending on stellar 
parameters. Negligible NLTE effects were obtained for Fe II. Fe I is subject to significant 
NLTE effects for low gravity (log g < 3) and very metal-poor models. An important 
consequence is that surface gravities of giants and very metal-poor stars derived by LTE 
analysis are in error with a magnitude strongly depending on log g/[Fe/H]. For example, 
LTE leads to a 0.26 dex lower gravity for T cff = 5000 K, log g = 2, and [Fe/H] = -3. 

4. Conclusions 

• Completeness of model atom for Fe I is important for a correct calculation of the 
Fe I/Fe II ionization equilibrium in the atmosphere of cool stars. 

• Thermalizing processes not involving electron collisions have to be included in the 
SE calculations for Fe I-II. Collisions with hydrogen atoms could be good candidates for 
such processes. 

• Fe I is affected by significant NLTE effects for giants and very metal-poor stars. 

• Only minor departures from LTE are obtained for Fe II. 
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